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Abstract. Exactly solvable models of ultracold Fermi gases are reviewed via their 
thermodynamic Bethe Ansatz solution. Analytical and numerical results are obtained 
for the thermodynamics and ground state properties of two- and three-component 
one-dimensional attractive fermions with population imbalance. New results for the 
universal finite temperature corrections are given for the two-component model. For 
the three-component model, numerical solution of the dressed energy equations confirm 
that the analytical expressions for the critical fields and the resulting phase diagrams at 
zero temperature are highly accurate in the strong coupling regime. The results provide 
a precise description of the quantum phases and universal thermodynamics which are 
applicable to experiments with cold fermionic atoms confined to one-dimensional tubes. 



PACS numbers: 02.30.1k, 03.75.Ss, 03.75.Hh, 64.70.Tg 
1. Introduction 

The experimental realization of Bose-Einstein condensates in dilute atomic gases [IH3] 
led to an explosion of ongoing research in ultracold matter physics. A remarkable 
development was the creation of the first molecular condensate in an ultracold degenerate 
Fermi gas (HE], after which the condensation of fermionic pairs was soon detected 
and shown to be a superfluid [6]. A relevant question in this fermionic context 
is if superfluidity can persist in a Fermi gas with imbalanced spin population. In 
principle, superfluidity may still occur in a mismatched case, and some theories with 
unusual pairings and exotic phases have been proposed, such as the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) phase, among others 018]. Subsequently the search for 
experimental confirmation of these new phases has been conducted by different groups 
[9j[T0]. Until now, however, only paired and polarized phases have been detected in 
three dimensions. 
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Further remarkable developments have involved the optical confinement of ultracold 
atoms to one dimension (ID), whereby atoms are trapped and cooled in an array of ID 
tubes. These systems include bosonic Rb atoms [THUS] an d fermionic 40 K atoms [13J. 
The most recent experimental breakthroughs involve ID Rb atoms in the attractive 
regime [H] and the realization of a ID spin-imbalanced attractive Fermi gas of 6 Li atoms 
under the degenerate temperature [15]. This experimental work in ID has highlighted 
the fundamental nature of exactly solvable models of quantum many-body systems. 
In particular, the Bethe Ansatz (BA) integrable models of Lieb and Liniger [16] for 
spinless bosons and of Yang [17] and Gaudin [18] for two-component fermions. The 
experimental work has further highlighted the deep and enduring significance of the 
BA [19J. More general BA integrable mult i- component fermions were initially studied by 
Sutherland [20J. In this paper we examine the BA integrable two- and three-component 
attractive ID Fermi gases with population imbalance. We shall see that these ultracold 
Fermi gases may be used to create nontrivial and exotic phases of matter. They also 
pave the way for the direct observation and further study of FFLO-like states. 

2. Two-component attractive Fermi gas with polarization 

2.1. The model 

We begin by reviewing the exactly solved two-component model, with Hamiltonian 



which describes N 5- interacting spin- 1 fermions of mass m constrained by periodic 
boundary conditions to a line of length L and subject to an external magnetic field H. 
The inter-component interaction g>iD can be tuned from strongly attractive to strongly 
repulsive via Feshbach resonance and optical confinement. The interaction is attractive 
for gn) < and repulsive for > 0. Here we will focus on the strongly attractive 
case, since this regime can be experimentally reached in ID [15] . For convenience, we 
define c = rugm/h 2 and a dimensionless interaction strength 7 = c/n for the physical 
analysis, with linear density n = N/L. 

The model (pQ), which exhibits SU(2) symmetry, was solved independently by 
Yang [T7] and Gaudin [18] using the nested BA. The energy eigenspectrum is given 
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where the quasimomenta {kj} of the fermions, satisfy the BA equations [T7J 

czn(ikL) TT ^- A ^ + ic / 2 
exp^L)- U 

i=i 3 1 1 

n X A a -k £ + i c/2 _ j-j A Q — + i c 
^ A a -A; £ -ic/2 = llA^-A^-ic' U 

for j = 1, . . . , N and a = 1, . . . , M, with M the number of spin-down fermions. Here 
{A a } are the rapidities for the internal spin degrees of freedom. 

The solutions to the BA equations ([3]) provide the ground state properties and 
the elementary excitations of the model. It was shown [2"Tjl2"2"] that the distribution 
of the quasimomenta in the complex plane for the ground state involves bound states, 
which can be interpreted as BCS-like pairs, and unpaired (excess) fermions. In the 
thermodynamic limit, i.e., N, L — > oo with N/L finite, and at zero temperature, all 
quasimomenta kj of N atoms form two-body bound states, i.e., kj = Aj ± i|c for 
j = 1, . . . ,M, accompanied by the real spin parameter Aj [21]. The BA equations ([3]) 
then become 

kjL =2-Klj + J2 e ( — j » j = 2M + 1,...,N (4) 



i=i 

N—2M / . i \ M 



j = l,...,M. 



where 6{x) = 2arctan2a; with Ij = —(N - 2M - l)/2, -(A - 2M - 3)/2, . . . , (N - 
2M- l)/2 and J 3 = —{M — l)/2, . . . , (M — 3)/2, (M— l)/2. 

Introducing the density of unpaired fermions p(fc) = dIj(k)/Ldk and the density of 
pairs o"(A) = dJj(A)/LdA it follows that 

1 1 y B |c|o-(A)dA 



2vr 2tt c 2 /4+ (jfc - A) 2 ' 
_ I l_ f B 2\c\a{A')dA' _ l_ r Q \c\p{k)dk 
CT( 7T ~ 2n J_ B c 2 + (A - A') 2 ~ 27 7_ Q c 2 /4 + (A - k) 2 (-)) 

in terms of which the total number of particles, the magnetization and the ground state 
energy per unit length are given by 

fB fQ 

n = 2 a(A)dA+ / p(k)dk, (6) 

J-B J-Q 

M z = ±(n t - ni ) = ± J p(k)dk, (7) 
Y= [ B (2A 2 - c 2 /2) a(A)dA + [ k 2 p(k)dk. (8) 

L J-B J-Q 

Here the integration boundary Q characterises the Fermi momentum in quasi- 
momentum space and B is the Fermi momentum in spin rapidity space. 
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Figure 1. Schematic representation showing the mixture of bound pairs and excess 
fermions in the strongly attractive regime. 



The BA equations (|5]) are essential to understanding the ground state properties 
of the model. The integration boundaries implicitly depend on the total number of 
particles and the magnetization. The strong coupling expansion for the ground state 
energy per unit length 
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(9) 



has been obtained by using the discrete BA equations (J3j) [22] or by iteration with the 
integral equations [23J. Here the polarization P is defined by P = (N^ — N±)/N and 
the binding energy is e& = Ft 2 n 2r ) 2 / '4m. The explicit relation between the polarization 
and the external field H is given further below in equation ( fl6l) . In this strong coupling 
region, the system behaves like a mixture of bound pair composites and excess single 
particles. The form of the ground state energy ([9]) reveals the nature of weak pair-pair 
attractive interaction and pair-unpair attractive interaction. A schematic representation 
of this mixture is depicted in Figure CD 



2.2. Thermodynamic Bethe Ansatz 

The thermodynamic properties as well as quantum phase transitions in this model 
system can be analyzed through the Thermodynamic Bethe Ansatz (TBA). This 
method was originally proposed by Yang and Yang [21] for spinless bosons and received 
important contributions for fermions by Takahashi [25], among others [26j . The basic 
idea of the method is that in the thermodynamic limit we can consider a continuous 
distribution function for the BA roots. The equilibrium state can be obtained by the 
condition of minimizing the Gibbs free energy G, 

G = E- HM Z -fxN- TS, (10) 

where /i is the chemical potential and S the entropy. 

This procedure gives rise to a set of coupled nonlinear equations, the TBA 
equations, for the dressed energies e b for paired and e u for unpaired fermions (see [25T - f2"T] 
for details) 

e b (k) = 2(k 2 -ix- -c 2 ) + Ta 2 * ln(l + e'^^) + Ta x * ln(l + e~^ T ), 
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e u (A;) = k 2 -fi--H + T ai * ln(l + e~^l T ) - T ^ a n * ln(l + ^{k)), 

TT OO 

In r] n (X) = ^ + a„ * ln(l + e~^ T ) + £ T nm * ln(l + ^(A)) (11) 

m=l 

where * denotes the convolution integral (/ * g)(X) = J™ f(X — X')g(X')dX' and the 
function 7/ n is the ratio of string densities, coming from Takahashi's hypothesis that 
complex quasimomenta group together to form strings in the complex plane [26]. The 
functions 

1 m\c\ 

«m A = — 12 

27T (m c/2) 2 + X 2 

and T mn (A) can also be found in Takahashi's book |26j . 

In the limit when T — > the TBA equations reduce to the dressed energy equations 

e b (A) = 2(A 2 - y, )-/ a 2 (A - A')e b (A / )rfA / - / fll (A - k)e u (k)dk, 



e u (k) = (k 2 - jj, ) - / ai (k - A)e h (A)dA 

2 J-B 



(13) 



The Gibbs free energy per unit length at zero temperature can be written in terms of 
the dressed energies as 



G(fi,H) = - [ e b {A)dA + — [ e u {k)dk, (14) 

TT J-B 27T J_q 

from which the density of fermions and the magnetization (or polarization P) follow via 
the relations 

- dG(ji, H)/dfi = n, -8G(p, H)/dH = m z = nP/2. (15) 

In general, the dressed energy equations provide a clear picture of band fillings with 
respect to the field H and the chemical potential \x at arbitrary temperatures [28J. 
These equations can be analytically solved in some special limits, such as in the strongly 
attractive regime |7| 3> 1 [27J discussed below. For comparison, we also present their 
numerical solution |29|. 



2.3. Phase diagram at zero temperature 



The set of equations ({TBI . (IT4"|) and (fT5l) were solved in [27, 29] for strongly attractive 
interaction using a lengthy iteration method. To leading order, the explicit form for the 
external field in terms of the density of fermions, the polarization and the interaction 
strength is given by 



h 2 n 2 
2m 
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— + 2P V 
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This equation determines the full phase diagram of the model and the critical values of 
the external field 



These results for if cl and H c2 are obtained from equation ( fl6|) by setting the polarization 
to P = and P = 1, respectively. The critical fields can be physically interpreted as 
separating three distinct quantum phases: (i) for H < H c i bound pairs populate the 
ground state, (ii) for H > H C 2 a completely ferromagnetic phase occurs, and (iii) for the 
intermediate range H c \ < H < H c2 paired and unpaired atoms coexist. 

Figure |2] illustrates the phase diagram in the n — H plane for the particular value 
|c| = 10. The dashed lines are plotted from equations ( TT7|) . The coloured phases are 
obtained by numerically integrating the dressed energy equations (fl3l) . The analytical 
results coincide well with the numerical boundaries. Here it is also clear that subject to 
the value of the external field, the system exhibits three quantum phases: a BCS-like 
fully paired phase; an unpaired, fully polarized phase; and a mixed, partially polarized 
phase, composed of BCS-like pairs and unpaired (excess) fermions. The mixed phase 
can be considered as a ID-analogue of the FFLO phase (see, e.g. [30H33] ). 

Calculations on the BA integrable model thus lead to a two-shell structure composed 
of the partially polarized FFLO-type phase in the centre of the trap surrounded by - 
depending on the strength of the external field - either fully paired or fully polarized 
wings [27J[3lJ[35]. This prediction was recently verified by the observation by Liao et 
al. [15] of three distinct phases in experimental measurements of ultracold 6 Li atoms 
in an array of ID tubes. A number of other authors have also investigated the BA 
integrable model in the context of the Fermi gas (see, e.g., [3&H39] ). 

2.4- TEA for low temperature 

In order to handle with the TBA equations ( fTTT) at low temperature, we can employ an 
expansion in terms of the polylogarithm function, which has so far been applied to the 
ID attractive Fermi gases of ultracold atoms up to order 1/c 2 (39HJT]. This approach 
is widely applicable to ID many-body systems with quadratic or linear bare dispersions 
in both the attractive and repulsive regimes. It is also expected that the analytical 
polylogarithmic function approach will play a central role in unifying the properties of 
attractive Fermi gases of ultracold atoms with higher symmetries. For example, the TBA 
equations for the ID Fermi gases with (^-function attractive interactions and internal 
spin degrees of freedom may be reformulated according to the charge bound states and 
spin strings characterizing spin fluctuations [ID] . Thus for strong attraction, the spin 




(17) 
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n 



Figure 2. Phase diagram for strong coupling |c| = 10 in the n — H plane (from 
Ref. [IS]). Good agreement is found between the analytical results (|17l) for the critical 
fields represented by the dashed lines and the numerical solution of the dressed energy 
equations (fT3|). given by the coloured phases. 



fluctuations that couple to non spin-neutral charge bound states are exponentially small 
and can be asymptotically calculated [40]. Thus the low energy physics is dominated 
by density fluctuations among the charge bound states. The full phase diagrams and 
thermodynamics of the ID attractive Fermi gases can be analytically calculated via this 
polylogarithm function approach. For spin- 1/2 attractive fermions, universal Tomonaga- 
Luttinger liquid behaviour was identified from the pressure given in terms of polylog 
functions (39]. Here we discuss some further results, on the universal nature of finite 
temperature corrections. 

For H 3> ksT, where ks is the Boltzmann constant, we can ignore the exponentially 
small corrections from the spin bound states and collect the terms up to I/7 2 in the 
TBA equations ( TTTT) . resulting in the simple form 

e b {k) « 2k 2 - A h {T,H), 

e u {k)^k 2 -A u {T,H), (18) 
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where 



A b (T, H) = 2 



2 c 



(p b (T,H)+4p»(T,H)) 



O 3 , (19) 



A U (T, JEO = Ai + Tln(2 cosh A) - ^p h (T, H) 



2T' 



1-B 4 



(20) 



Here we have denoted B = taring and J = -^p u (T, H). In order to examine the low 
temperature Luttinger liquid signature, we can safely ignore the exponentially small 
contribution from the spin wave bound states in the function A u . Thus integration by 
parts gives the effective pressures 

y/2 f°° Jldt 



p\T,H) 



P U (T,H) 



w 2 h 2 

2m 



A b (T,H) ' 
k B T 



K 2 h 2 
2m 



1 + e 
y/ede 

-A U (T,H) 



(21) 



1 + e" 



The integrals in (I2T1) can be calculated explicitly using Sommerfeld expansion at low 
temperatures. We assume that there exist two Fermi 1.6. , db Fermi sea of bound 

pairs with a cut-off potential A h (T, H)/2 and a Fermi sea of unpaired fermions with 
a cut-off potential A n (T, H). After some lengthy iteration with the relations ( fl5l) we 
obtain the free energy 



F(T,H) 



h 2 ir 2 n 3 P 




-nPH 
2 



nil - P) — , 

2m K 4 



(22) 



where /iJJ 



h 2 -w 2 n 2 P 2 



and /^q 



h 2 7T 2 n 2 {l-P) 2 



2m 32m 

to the ground state energy ([9]) as T — > 0. 



. We see clearly that the free energy reduces 



For the temperature k^T <C Ep, where Ep is the Fermi energy, the spin wave 
fluctuation is frozen out due to the large magnetic field. The leading low temperature 
correction to the free energy (122]) is 



F(T,H) = E (H) 



7iCk 2 B T 2 



6H 



1 1 



(23) 
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which belongs to the universality class of the Gaussian model with central charge 
C = 1 [42J. In the above equation, the ground state energy E (H) is as given in 
(EJ), subject to an additional term —nPH/2. The group velocities for bound pairs and 
unpaired fermions are 



respectively Here the Fermi velocity is Vp = hirn/m. The above result indicates 
that the low energy physics for ID strongly attractive fermions in the gapless phase 
can be described by a two-component Tomonaga-Luttinger liquid model as long as the 
ferromagnetic spin-spin interaction is frozen out. 

In order to recognize a Tomonaga-Luttinger liquid signature in the gapless phase 
in the low temperature limit, we examine the temperature dependent relations for the 
magnetization. In general [43J, one should expect a magnetization minimum due to a 
crossover from a Tomonaga-Luttinger liquid with a linear dispersion to a state governed 
by the nonrelativistic dispersion e oc k 2 . This magnetization minimum does exist in 
the gapless phase in attractively interacting fermions. We plot the magnetization vs 
temperature from the free energy (122"]) in Figure El where we observe a clear minimum 
of the magnetization for different magnetic fields. This comes about due to a crossover 
from the hardcore bosonic signature of the Tomonaga-Luttinger liquid to polarized free 
fermions. For further understanding this signature, we derive the magnetization from 
( )23l) where we consider the low temperature limit, i.e., T — > (in natural units) 



low temperatures depends on the signs from the term in the brackets of equation (T25]) . 
This part indicates the existence of a minimum of the magnetization. In the gapless 
phase two Fermi liquids are coupled through pair-unpaired fermion scattering. The 
linear field dependent magnetization is a consequence of the fact that the total number 
of fermions is fixed. 

So far we have considered the simplest ID exactly solvable model in the scenario of 
ultracold Fermi gases. Generalizations to three and more components can be performed 
and are discussed in the next sections. 




(24) 




(25) 



where = dE Q (H) /OH . The finite temperature contribution to the magnetization at 
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Figure 3. Magnetization m z vs temperature at different externai field values (in the 
units 2m = H = 1). The minimum of the magnetization is clearly observed for each 
field value. 



3. Three-component attractive Fermi gas with polarization 



3.1. The model 



To describe a three-component Fermi gas, the same type of Hamiltonian ([I]), with a 
kinetic and a contact potential interacting term, can be considered. However, now the 
N fermions can occupy three possible hyperfine levels (|1), |2) and |3)) and the Zeeman 
term is expressed in terms of two external fields. The Hamiltonian reads [20] 

F =-£EJ^ + ^ £ ^-x^ + Y.N'eWs.B). (26) 

i=l * l<i<j<N i=l 

The last term denotes the Zeeman energy, where N % is the number of fermions in state 
\i) with Zeeman energy e % z determined by the magnetic moments [i B and the magnetic 
field B. This term can also be written as —Hi(N l - N 2 ) - H 2 {N 2 - N 3 ) + Ne, where 
the unequally spaced Zeeman splitting in three hyperfine levels can be specified by two 
independent parameters Hi = e — €z(jj, b , B) and H 2 = €%{^ B , B) — e, with e the average 
Zeeman energy. We use the same notation and conventions as in the previous case. 

The Hamiltonian (|26l) exhibits SU(3) symmetry and was solved by Sutherland in 
the sixties by means of the nested BA [20] . The energy eigenspectrum is again given in 
terms of the quasimomenta {k{} of the fermions by (T5]), but now satisfying [201121] 

a t\ TT%-A^ + ic/2 
e Xp (i%L) = n 

i=i 3 1 1 

n A a - h + i c/2 tt A a — 4- 1 c -p- r Aq, — — i c/2 
, a A a -h-ic/2~ L} x A a - A p - ic\\A a - X e + ic/2' 
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n' h y- Ai + ic/2 = -tt A m - A^ + ic 
;=i A,-A,-ic/2 l = lA M -A,-ic 1 J 

for j = 1, . . . , N, a — 1, . . . , Mi, n — 1, . . . , M 2 , with quantum numbers Mi = A^ + 2A" 3 
and M 2 = iV 3 . The parameters {A a , A m } are the rapidities for the internal hyperfme 
spin degrees of freedom. For the irreducible representation [3 Ns 2 N2 l Nl \ , a three-column 
Young tableau encodes the numbers of unpaired fermions (Aq = A^ 1 — A^ 2 ), bound pairs 
(A" 2 = A^ 2 — A^ 3 ) and trions (A" 3 = N 3 ). 

The ground state energy per unit length 

E ir 2 n 3 / 8n 2 + 4n 3 \ n 2 c 2 



1 + 

7i 2 nl ( 12m + 6n 2 + 16n 3 \ 2 

+ — + 3M )~ 2n3C 

+ fr4^ 1+ 12m + 32n 2 + 18n 3 ^ (2g) 

was obtained in [H] by solving the BA equations (|27|) . Here n a = N a /L, a = 1,2,3 
are the densities of unpaired fermions, bound pairs and trions, respectively, with the 
constraint n = m + 2n 2 + 3n 3 . 

The phase diagram of the system for the strongly attractive interaction regime can 
also be obtained using the dressed energy formalism. For this model, a richer scenario, 
with more quantum phases, is expected. Indeed, as we shall discuss further below, a 
trion phase, which consists of three-body bound states, along with a number of mixed 
phases, are found 



3.2. Dressed energy formalism and phase diagrams at zero temperature 

In the thermodynamic limit the TBA equations are found by minimizing the Gibbs free 
energy G = E + n\H\ + n 2 H 2 — fiN — TS. In the limit T — > 0, the dressed energy 
equations obtained [261IHII45] are 

e ( 3 )(A) = 3A 2 - 2c 2 - 3/i - a 2 * e (1) (A) 

-[a 1 + a 3 ]*e (2) (A)-[a 2 + a 4 ]*e (3) (A) 

e (2) (A) = 2A 2 - 2/x - °- - H 2 - m * e (1) (A) 

-a 2 *e 2 (A)-[ai + a 3 ]*e (3) (A) (29) 
e W(A;) =k 2 -fi-H 1 -a 1 * e {2 \k) - a 2 * e {3) (£;). 

Here are the dressed energies and aj*e^(x) = J^q" CLj(x — y)e^ a \y)dy. The negative 
part of the dressed energies e^(x) for a; < \Q a \ correspond to the occupied states in the 
Fermi seas of trions, bound pairs and unpaired fermions, with the positive part of 
corresponding to the unoccupied states. The integration boundaries Q a characterize the 
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"Fermi surfaces" at e^ a \±Q a ) = 0. The zero-temperature Gibbs free energy per unit 
length can be written in terms of the dressed energies as 



3 n r +Qa 

G = J2~ e^{x)dx, (30) 

a=l 77 J-Qa 

from which physical quantities are obtained through the thermodynamic relations 
dG dG dG 

The dressed energy equations (1291) can be analytically solved just in some special 

limits, such as the strongly attractive coupling regime. This was discussed in [44] . where 

the analytical expressions 

2 2 f 4m 8n 2 AnA IOtt 2 ^ 

Hi = 7T TI-, 1 — + — r + -r-r H r-r - 

1 V 9|c| |c| |c| J 27\c\ 

■K 2 n\ ( 4m 32n 2 4n 3 \ 2c 2 
1 H — ; — r H r-r + ' 



3|c| 9|c| 3|c|/ 3 
7r 2 no / 4m 40n 2 16n s \ 167r 2 n 3 



27|c| 3|c| / 9|c| 
2vr 2 n 2 / 4m 32n 2 8n 3 \ 5c 2 



9 \ 3|c| 9|c| 9|c|/ 6 
for the energy transfer relations were found (in units of h 2 /2m). These equations 
determine the full phase diagram and the critical field values activated by the fields 
Hi and H 2 . 

Here we numerically solve these equations to confirm the analytical expressions for 
the physical quantities and the resulting phase diagrams of the model. Basically, for each 
value of the integration boundaries Q a , a = 1,2,3, the equations (12"9"1) are converted into 
a finite size single matrix equation that can be solved using the condition e^(±Q a ) = 
to get the dressed energies, the fields and the chemical potential. Then the Gibbs free 
energy is obtained and from it the polarizations and the linear density are found. Each 
phase can be identified by properly separating the input in distinct sets containing 
different combinations of vanishing/ non vanishing integration boundaries. For example, 
in the pure unpaired phase we have {Q 2 = Q3 = 0,Qi 7^ 0}, in the mixed phased 
composed of unpaired fermions and bound pairs {Q3 = 0, Qi 7^ 0, Q2 7^ 0} and so forth. 
There are seven different combinations and consequently, seven different phases, which 
are computed and the results collected to generate the complete phase diagrams. 

Figures H] and [5] show the polarizations n\jn and n 2 /n in terms of the fields Hi 
and H 2 . There are three pure phases: an unpaired phase A, a pairing phase B and a 
trion phase C and four different mixtures of these states. A good agreement is observed 
between the analytical predictions obtained from equations ( 1321) represented by black 
lines and the numerical solutions obtained by integrating the dressed energy equations 
(129]) represented by white dots. 
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Figure 4. Phase diagram showing the polarization ri\/n versus the fields Hi and 
H% for strong interaction with \c\ = 10 and n = 1. There are three pure phases: an 
unpaired phase A, a pairing phase B and a trion phase C and four different mixtures 
of these states. The black lines are plotted from the analytical results (1521) while the 
white dots correspond to the numerical solutions of the dressed energy equations (|29[) . 
The numerical phase transition boundaries coincide well with the analytical results. 




77.5 



Figure 5. Phase diagram showing the polarization 77,2 /n versus the fields Hi and 
H 2 for strong interaction with |c| = 10 and n = 1. The black lines plotted from the 
analytical results p2|) are in good agreement with the numerical solutions (white dots) 
of the dressed energy equations (|2^|) . 



The ground state energy versus the fields Hi and H 2 can be determined from the 
ground state energy ( )28l) with the densities n\ and n 2 obtained from equation (I32j) . 



Exactly solvable models and ultracold Fermi gases 



14 




Figure 6. Ground state energy vs Zeeman splitting for strong interaction \c\ = 10 
and n = 1. Good agreement is found between the analytical results (I32[) represented 
by black lines and the numerical solutions obtained by integrating the dressed energy 
equations (1^1) . represented by white dots. 



Figure [6] shows the energy surface for all possible phases shown in Figures H] and |5j 
This figure demonstrates the interplay between different physical ground states. For 
certain values of Hi and H^, a mixture of unpaired fermions, BCS-like pairs and trions 
(A + B + C) populates the ground state. 

For low temperature, a similar investigation employing the polylog function can 
also be performed |41j . 

4. Conclusion and perspectives 

We have examined the two- and three-component attractive ID Fermi gases with 
population imbalance via their TBA solution. For the two-component model, we 
reviewed the strong coupling expansion and the identification of quantum phases. New 
results for the universal finite temperature corrections were also discussed. For the 
three-component model, numerical solution of the dressed energy equations confirm 
that the analytical expressions for the critical fields and the resulting phase diagrams 
at zero temperature are highly accurate in the strong coupling regime. Both models 
exhibit rich phase diagrams with a variety of quantum phases. Just as the three distinct 
phases of the two-component model - the BCS-like paired, fully polarized and partially 
polarized phases - have been detected by Liao et al. [15] in a recent experiment with 
ultracold 6 Li atoms in an array of ID tubes, it is to be hoped that the more exotic 
phases of the three-component model - including the trion phase - will be detectable in 
future experiments. 
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From the experimental point of view [TB], the array of ID tubes is created within 
2D optical lattices. In order to make the lowest transverse mode populated in each tube, 
the thermal energy /c^T and the Fermi energy are required to be much smaller than 
the transverse confinement energy. In this sense, the system is well controllable and 
practicable only for the strongly attractive regime in the quasi-lD trapping. Precisely 
in this regime, the analytical results may provide direct application to fitting the 
experimental data, such as the density profiles and phase diagram |16] . 

Mult i- component Fermi gases with more than three species can also be trapped 
and manipulated. For this type of Fermi gas, bound multi-body clusters are expected 
to appear above certain critical interaction strengths |I7]. The thermodynamic 
properties and phase diagrams of ID attractive multicomponent Fermi gases can also 
be investigated through the solvable models exhibiting SU(N) symmetry [ID]. Closed 
form expressions for the thermodynamics and equation of state of such models will 
provide further insight towards understanding the nature of many-body effects and 
different pairing states with higher spin symmetry. The trapping potentials can be 
accommodated into the equation of state within the local density approximation. It 
is clear that the Bethe Ansatz will continue to prosper as an essential tool for their 
description. 
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